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Abstract  In  1884,  the  International  Meridian  Conference 
recommended  that  the  prime  meridian  “to  be  employed  as  a 
common  zero  of  longitude  and  standard  of  time-reckoning 
throughout  the  globe”  pass  through  the  “centre  of  the  transit 
instrument  at  the  Observatory  of  Greenwich”.  Today,  tourists 
visiting  its  meridian  line  must  walk  east  approximately  102  m 
before  their  satellite-navigation  receivers  indicate  zero  lon¬ 
gitude.  This  offset  can  be  accounted  for  by  the  difference 
between  astronomical  and  geodetic  coordinates — deflection 
of  the  vertical — in  the  east-west  direction  at  Greenwich,  and 
the  imposed  condition  of  continuity  in  astronomical  time. 
The  coordinates  of  satellite-navigation  receivers  are  pro¬ 
vided  in  reference  frames  that  are  related  to  the  geocentric 
reference  frame  introduced  by  the  Bureau  International  de 
I’Heure  (BIH)  in  1984.  This  BIH  Terrestrial  System  provided 
the  basis  for  orientation  of  subsequent  geocentric  reference 
frames,  including  all  realizations  of  the  World  Geodetic  Sys¬ 
tem  1984  and  the  International  Terrestrial  Reference  Frame. 
Despite  the  lateral  offset  of  the  original  and  current  zero- 
longitude  lines  at  Greenwich,  the  orientation  of  the  meridian 
plane  used  to  measure  Universal  Time  has  remained  essen¬ 
tially  unchanged. 
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1  The  situation  at  Greenwich 

The  Airy  Transit  Circle  is  a  nineteenth-century  telescopic 
instrument  at  the  Royal  Observatory,  Greenwich,  famous 
as  the  origin  for  global  longitude.  Recognized  as  a  place 
of  world  heritage  by  UNESCO,  the  observatory  is  now  a 
museum  with  a  line  in  the  pavement  running  north  to  south, 
and  with  signage  marking  it  as  the  “Prime  Meridian  of  the 
World”.  Today,  its  longitude  is  00°00'05.3"W  (Howse  1997) 
in  the  International  Terrestrial  Reference  Frame  (ITRF),  and 
in  the  World  Geodetic  System  (WGS  84)  used  by  the  U.S. 
Global  Positioning  System  (GPS).  The  observatory’s  non¬ 
zero  longitude  is  widely  noticed  by  tourists  with  satellite- 
navigation  receivers,  and  by  people  using  software  to  display 
geo-referenced  imagery  or  maps  of  the  Earth.  The  ground 
trace  of  the  WGS  84  and  the  ITRF  zero-longitude  meridian 
plane  is  located  approximately  102  m  east  of  the  telescope 
(Fig.  1).  The  apparent  discrepancy  raises  questions  of  how 
and  when  it  arose,  and  whether  a  worldwide  system  of  lon¬ 
gitude  has  been  systematically  shifted  by  this  same  amount. 

In  Sects.  2  and  3,  the  authors  show  that  the  deflection  of 
the  vertical  (DoV)  can  account  for  the  entire  longitude  shift 
at  Greenwich.  In  Sect.  4,  the  authors  also  compare  astrogeo- 
detic  and  gravimetric  DoVs  and  their  uncertainties,  using  a 
recent  global  gravitational  model,  for  a  number  of  astronom¬ 
ical  observatories  that  contributed  to  past  determinations  of 
UT1 .  This  provides  information  about  the  overall  orientation 
between  the  former  (astronomical)  and  modern  (geodetic) 
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Fig.  1  The  Airy  meridian  ( dotted  line )  and  the  ITRF  zero  meridian  {solid  line).  Imagery©  2014  Google  Maps,  Infoterra  Ltd  &  Bluesky 


longitude  systems.  Because  books  for  the  general  public  (e.g., 
Malin  and  Stott  1989;  Howse  1997;  Jennings  1999;  Dolan 
2003;  Murdin  2009),  as  well  as  the  web  sites  of  the  UK’s 
National  Maritime  Museum  (Sinclair  1997)  and  Wikipedia 
(IERS  Reference  Meridian  2014),  either  ignore  the  longitude 
offset  at  Greenwich,  acknowledge  it  but  do  not  assign  a  cause, 
or  provide  incorrect  accounts,  the  authors  also  consider  some 
of  the  conjectured  explanations  in  Sect.  5. 


2  Methods  of  determining  longitude  and  Earth 
rotation 

Traditional  optical  instruments  for  measuring  Earth  rotation 
included  visual  and  photographic  zenith  tubes,  circum- 
zenithals,  Danjon  astrolabes,  and  transit  telescopes,  such  as 
the  Airy  Transit  Circle.  These  instruments  measured  local 
sidereal  time  by  the  passage  of  specially  designated  “clock 
stars”  across  a  plane  of  reference  established  in  part  by  the 
local  direction  of  gravity.  These  measurements  were  realized 
with  respect  to  the  local  vertical  using  a  basin  of  mercury,  and 
the  astronomical  latitudes  and  longitudes  of  the  observatories 
were  thus  affected  by  local  gravity,  because  their  coordinates 
were  determined  from  the  same  observations. 

The  sidereal  time  observed  from  these  instruments  was 
converted  to  mean  solar  time  according  to  a  conventional 
relationship,  and  time  signals  were  adjusted  to  match  the 
astronomical  observations.  The  difference  between  astro¬ 
nomically  determined  local  mean  times  t  (in  hours)  is  pro¬ 
portional  to  the  difference  in  local  longitudes  A  (in  degrees): 

t-t0  =  (A-A0)/15,  (1) 

where  Aq  and  to  refer  to  the  astronomical  longitude  and  time 
of  the  reference,  or  “prime”,  meridian.  The  1884  Interna¬ 


tional  Meridian  Conference  recommended  that  the  prime 
meridian,  Ao  =  0,  pass  through  the  Airy  Transit  Circle 
at  Greenwich  (Explanatory  Supplement  1961),  establish¬ 
ing  the  mean  solar  time  determined  at  the  Airy  Transit 
Circle — Greenwich  Mean  Time  (GMT) — as  an  internation¬ 
ally  recognized  basis  for  global  timekeeping  and  navigation. 
Any  astronomical  time  service  could  thereby  convert  its  own 
determination  of  local  mean  time  to  GMT  by  removing  its 
longitude  difference  from  Greenwich. 

Because  of  significant  discrepancies  between  the  earli¬ 
est  radio  time-signal  broadcasts,  the  Bureau  International 
de  VHeure  (BIH)  was  established  at  the  Paris  Observatory 
in  1912  to  help  coordinate  the  national  time  bureaus  and 
improve  longitudes  globally  (Guinot  2000).  For  purposes  of 
global  timekeeping  by  the  BIH,  the  Greenwich  Observatory 
was  superseded  in  1929  by  a  statistical  construct  known  as 
the  “mean  observatory”,  of  which  Greenwich  became  one  of 
many  contributors.  The  formation  of  GMT,  now  called  Uni¬ 
versal  Time,  from  this  “mean  observatory”  involved  adopted 
corrections  to  the  contributed  data  and  weights  that  could 
change  annually  (Feissel  1980).  To  maintain  a  continuous  and 
consistent  series  of  Universal  Time,  the  BIH  also  assumed 
the  responsibility  for  establishing  and  refining  the  terrestrial 
reference  systems  used  for  the  analyses  (Guinot  2000).  The 
early  BIH  terrestrial  reference  frames  were  entirely  based  on 
astronomical  coordinates  (Mueller  1969). 

Delays  in  the  analysis  of  polar  motion  from  the  Inter¬ 
national  Latitude  Service  (ILS)  made  it  impractical  for  the 
BIH  to  correct  Universal  Time  for  the  motion  of  the  rota¬ 
tional  pole  until  the  middle  of  the  twentieth  century  (Guinot 
2000).  Universal  Time  observed  from  a  given  observatory 
was  designated  UT0,  whereas  the  distributed  astronomical 
time  scale  corrected  for  polar  motion  became  known  as  UT1 . 
After  the  BIH  adopted  various  reference  poles  of  different 
epochs  during  the  1950s  and  1960s  (Robbins  1967),  the  BIH 
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standardized  on  the  average  ILS  pole  from  1900-1905 — 
known  as  the  Conventional  International  Origin — as  part 
of  an  update  to  celestial  and  terrestrial  reference  systems. 
To  preserve  the  continuity  of  UT1-UT0  for  contributing 
stations  within  this  new  “1968  BIH  System”,  the  adopted 
longitudes  of  observatories  were  systematically  adjusted  for 
the  change  of  pole.  Although  the  Airy  Transit  Circle  was 
no  longer  operational  as  a  BIH  station,  its  astronomical  lon¬ 
gitude  would  have  changed  from  zero  to  0.2927"W  (about 
6  m)  within  this  new  system  (MacDonald  1985;  Schmid 
1974). 

Beginning  in  1973,  the  BIH  began  incorporating  estimates 
of  polar  motion  from  the  U.S.  Navy  TRANSIT  (satellite 
Doppler)  navigation  system  into  its  own  processes  (Guinot 
1979).  Other  experimental  techniques,  such  as  Very  Long 
Baseline  Interferometry  (VLB I),  connected-element  radio 
interferometry,  and  lunar  laser  ranging  (LLR)  were  gradually 
introduced  along  with  the  traditional  optical  measurements 
to  estimate  UT1.  By  1984,  the  BIH  discounted  optical 
astrometry  completely,  and  established  a  new  terrestrial  coor¬ 
dinate  system  designated  as  the  “BIH  Terrestrial  System”,  or 
BTS  84  (BIH  1985).  This  system  was  based  on  weighted 
solutions  of  Earth-orientation  parameters  from  analysis  cen¬ 
ters  and  networks  of  VLBI,  SLR  and  LLR  geodetic  stations 
that  were  independent  of  the  former  network  of  optical  sys¬ 
tems  (Boucher  and  Altamimi  1985,  1986).  In  constructing 
BTS  84,  the  solution  of  network  orientation  (and  thus  its 
corresponding  meridian  of  zero  longitude)  was  constrained 
to  match  the  ensemble  of  Earth-orientation  measurements 
made  from  1980  through  1983  relative  to  atomic  time¬ 
keeping.  Thus,  the  series  of  UT1  determinations  based  on 
non-traditional  techniques  maintained  continuity  with  the 
former  series  dominated  by  optical  astrometry  (within  the 
uncertainty  of  the  optical  measurements). 

By  1988,  the  BIH  and  the  International  Polar  Motion 
Service  (successor  to  the  ILS)  were  both  replaced  by  the 
International  Earth  Rotation  Service  (IERS)  (Wilkins  2000), 
which  assumed  responsibility  for  maintaining  the  terrestrial 
and  celestial  reference  frames  consistent  with  the  Earth- 
orientation  parameters  it  published  routinely.  The  subsequent 
terrestrial  reference  frames  of  the  IERS — in  particular,  the 
various  realizations  of  the  International  Terrestrial  Refer¬ 
ence  Frame  (ITRF) — preserve  the  orientations  relative  to  the 
reference  meridian  of  the  BTS  (see  Jekeli  2006,  pp.  3-14). 
Modern  geodetic  reference  frames,  such  as  the  ITRF,  are 
defined  by  the  Cartesian  position  and  velocity  coordinates 
of  an  ensemble  of  stations  in  three  dimensions.  Within  this 
context,  the  corresponding  reference  (zero)  meridian  is  the 
Earth-fixed  X-Z  plane,  where  longitude  is  measured  around 
the  Z-axis  eastward  from  that  plane.  The  IERS  Conventions 
(IERS  2010)  recommend  the  term  ‘ITRF  zero  meridian’  for 
this  plane,  although  the  terms  ‘IERS  Reference  Meridian’ 


and  ‘International  Reference  Meridian’  (IRM)  are  also  in 
use. 


3  Astrogeodetic  and  gravimetric  DoV 
determination 


Natural  coordinates  [0,  A,  H]  are  related  to  geodetic  lati¬ 
tude,  ip,  longitude,  k,  and  ellipsoidal  height  h ,  by  (Heiskanen 
and  Moritz  1967): 


£  =  &  ~  V  (a) 
rj  =  (A  —  k)  cos  ip  (b) 
N  —  h  —  H  (c) 


(2) 


where,  §  is  the  meridional  component  of  the  DoV,  rj  the  prime 
vertical  component  of  the  DoV  and  N  is  the  geoid  undula¬ 
tion.  These  equations  relate  the  geodetic  quantities  [ip,  k,h] 
to  the  physical  quantities  [0 ,  A,  H]  associated  with  the  grav¬ 
ity  field.  If  the  geodetic  coordinates  [ip,  k,  h\  of  a  point  have 
been  determined  with  respect  to  a  given  datum,  and  if  the 
natural  coordinates  [0,  A,  H]  have  been  determined  also 
(e.g.,  from  astronomical  observations  and  from  spirit  level¬ 
ing),  then  one  can  obtain  the  deflection  components  (Castro, 
Castro)  and  the  geoid  undulation  ( N )  at  that  point  from  Eq. 
(2).  Today,  such  a  geometric  determination  can  be  realized 
by  making  astronomical  observations  of  0  and  A  at  a  point 
whose  geodetic  coordinates  [ip,  k,  h]  have  been  determined 
from  GNSS  positioning,  and  whose  orthometric  height,  H, 
has  been  determined  from  spirit  leveling.  Such  a  geometric 
determination  yields  astrogeodetic  DoV  components  that  are 
specific  to  the  particular  datum  to  which  the  geodetic  coor¬ 
dinates,  ip  and  k,  refer. 

The  DoV  components  (§,  q)  also  represent  the  slopes  of 
the  geoid  surface  with  respect  to  the  surface  of  an  equipo- 
tential  reference  ellipsoid  along  the  meridian  and  the  prime 
vertical,  respectively,  (Heiskanen  and  Moritz  1967,  section 
2-22): 

§grav  =  -^  (a) 

Vv  =  -^alWj’ 

where  R  is  a  mean  Earth  radius.  The  sign  convention  in  Eqs. 
(2)  and  (3)  is  such  that  §  is  positive  when  the  astronomi¬ 
cal  zenith,  Za,  is  north  of  the  geodetic  zenith,  Zg;  and  r\  is 
positive  when  Za  is  east  of  Zg  (Torge  2001,  p.  219).  These 
slopes  can  be  determined  gravimetrically,  either  using  the 
integral  formulas  of  Vening  Meinesz,  or  using  a  spherical 
harmonic  representation  of  the  Earth’s  gravitational  poten¬ 
tial.  The  gravimetric  determination  of  DoV  requires  dense 
measurements  of  gravity  and  detailed  mapping  of  the  topog¬ 
raphy  near  the  computation  point  for  the  effective  use  of  the 
integral  formulas,  or,  equivalently,  a  very  high  degree  and 
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Fig.  2  The  reference  ellipsoid  from  a  point  on  its  minor  axis  above 
the  north  pole.  O  =  spheroid  axis  passing  through  the  geocenter,  Eq 
=  ground  level  equipotential  surface,  Sph  =  spheroidal  parallels,  OA'= 
geodetic  meridian  plane  at  K,  SA  =  astronomical  meridian  plane  at  K 
(Bomford  1980) 


order  gravitational  field  model.  EGM2008  is  a  gravitational 
model  that  is  complete  to  spherical  harmonic  degree  2190  and 
order  2159  (Pavlis  et  al.  2012)  and  carries  enough  resolving 
power  to  support  the  determination  of  gravimetric  deflec¬ 
tions  with  adequate  accuracy,  to  estimate  the  astronomical 
coordinates  for  the  Airy  Transit  Circle  by: 

^grav  =  <P  +  §grav  (a)  ^ 

Agrav  =  A  ?7grav  '  Sec  (^) 

These  “synthetic”  astronomical  coordinates  can  then  be 
compared  to  independent  astronomical  observations.  If  the 
orientation  of  the  zero  meridian  plane  is  the  same  for  both 
astronomical  and  geodetic  longitudes,  Eqs.  (2b)  and  (3b)  will 
yield  consistent  values  for  r\  (Bomford  1980,  p.  100). 

The  authors  used  the  EGM2008  spherical  harmonic  coef¬ 
ficients  to  degree  2190,  and  calculated  the  deflections  at 
Greenwich  via  harmonic  synthesis,  to  obtain  §grav  =  2.156" 
and  r/grav •  sec  ip  =  5.502"  assuming  zero  elevation,  a  result 
which  agrees  with  the  estimates  of  2.15"  and  5.51"  made  by 
Ekman  and  Agren  (2010)  via  numerical  differentiation  of  the 
gridded  version  of  the  EGM2008  geoid.  The  results  predict 
an  astronomical  longitude  of  Agrav  =  00°  00^0. 19"E  for  the 
Airy  Transit  Circle.  Within  the  estimated  EGM2008  com¬ 
mission  error  in  77grav  sec  ip ,  zb  0.47"  1  a  (see  also  Pavlis  et  al. 
2012,  section  5),  the  predicted  astronomical  longitude  is  in 
good  agreement  with  its  originally  adopted  value  of  zero. 

Figure  2  illustrates  the  reference  ellipsoid  as  seen  from 
a  vantage  point  along  its  minor  axis  above  the  north  pole 


(point  O).  The  line-segment  QC  identifies  the  plane  of  the 
astronomical  reference  meridian  of  Greenwich,  A  =  0.  It  is 
parallel  to  O  C",  the  plane  of  the  geodetic  reference  meridian 
of  G",  A  =  0.  The  angle  C'OC "  is  exactly  77-  sec  ip ,  the  east- 
west  component  of  the  DoV  at  Greenwich,  projected  onto 
the  equator.  As  illustrated,  the  shift  eastward  at  Greenwich 
represents  only  a  lateral  transfer  of  the  trace  of  the  geodetic 
longitude  A  =  0  meridian  across  the  surface  of  the  Earth. 
Thus,  the  geodetic  prime  meridian  at  Greenwich  has  the  same 
orientation  as  a  function  of  time  as  the  astronomical  meridian. 
In  order  for  the  plane  of  the  geodetic  prime  meridian  to  pass 
through  the  geocenter,  its  trace  on  the  ground  in  the  vicinity 
of  Greenwich  must  move  to  the  east  by  about  102  m,  the  sign 
and  magnitude  of  which  is  correctly  predicted  by  EGM2008. 


4  Astrogeodetic  and  gravimetric  DoV  comparisons 

The  astrogeodetic  deflection  ^astro  from  Eq.  (2b)  can  be 
expected  to  agree  with  its  corresponding  gravimetric  value 
r\ grav  only  if  the  astronomical  and  geodetic  origins  of  longi¬ 
tude  are  parallel  planes.  This  condition  is  not  guaranteed  in 
principle  (see  also  Tscherning  1986),  but  it  can  be  tested  at 
surface  points  where  gravimetric  determinations  of  77  have 
been  made,  and  where  astronomical  and  geodetic  longitudes 
A  and  A  are  both  determined.  The  difference  A  77  =  77astro  - 
77 grav  includes  errors  in  the  determination  of  the  coordinates, 
as  well  as  errors  of  commission  and  omission  in  the  gravi¬ 
metric  deflection.  For  a  meaningful  test,  these  errors  should 
be  small  relative  to  the  magnitude  of  any  supposed  misalign¬ 
ment  of  the  longitude  origin. 

A  slightly  different,  but  equivalent,  approach  follows  from 
Ekman  and  Agren  (2010),  who  “synthesized”  astronomical 
coordinates  per  Eq.  (4),  and  then  compared  values  of  Agrav 
to  the  historical  astronomical  longitudes  (A).  For  fundamen¬ 
tal  observatories  at  Stockholm,  Kpbenhavn,  and  Greenwich, 
their  results  indicated  differences  between  A  and  Agrav  of 
0.4",  0.9",  and  0.2",  respectively.  Here,  the  authors  extended 
this  approach  by  correcting  the  conventional  BIH  coordinates 
of  astronomical  stations  before  1984  (1968  BIH  System)  for 
local  deflection,  and  then  comparing  the  corrected  coordi¬ 
nates  to  their  geodetic  coordinates  with  respect  to  WGS  84. 
If  the  origin  for  longitude  was  the  same  both  before  and  after 
1984,  then  the  longitude  residual  A  -  Agrav  for  each  BIH 
optical  station  should  be  small  relative  to  the  uncertainty  of 
its  estimated  deflection,  cr rj -sec  y- 

Consideration  was  limited  to  optical  stations  that  could  be 
confidently  identified,  and  also  contributed  to  the  BIH  deter¬ 
mination  of  UT1  prior  to  BTS  84  (Table  1).  Where  geodetic 
coordinates  were  unavailable,  the  authors  resorted  to  satel¬ 
lite  imagery  available  via  virtual  globe  services  (e.g.,  Google 
Earth™,  Bing™  Maps,  etc.)  to  approximate  the  geodetic  lon¬ 
gitude  of  extant  shelters  or  their  ruins.  For  stations  that  were 


Springer 


Table  1  EGM2008 -implied  versus  Astrogeodetic  longitude  differences  for  BIH-contributing  optical  stations 
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72-83 

70-83  Pecny  49°54'50.5"N  14°47'12.1"E  549  8.4  8.1  0.90 

68-78  Richmond  25°36'50.8"N  80°23'03.0"W  8  7.2  7.8  0.42 
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co-located  with  geodetic  observatories,  the  precise  ITRF 
coordinates  of  nearby  GPS  or  SLR  stations  confirmed  that 
geolocation  errors  of  the  virtual  globes  were  at  the  level  of  a 
few  meters,  which  is  generally  below  the  uncertainty  in  /7grav 
for  EGM2008.  In  most  cases,  affiliated  personnel  provided 
corroboration  of  extant  structures,  or  geodetic  coordinate 
estimates  based  on  historic  local  surveys  or  GPS  readings; 
this  was  the  only  means  of  positioning  some  long-demolished 
stations. 

For  the  analysis,  astronomical  longitudes  were  adopted 
from  BIH  annual  reports  between  1968  and  1983,  except  that 
of  Washington  (W),  which  was  adopted  from  Hughes  et  al. 
(1975).  Small  annual  and  semi-annual  coordinate  corrections 
estimated  by  BIH,  which  could  have  been  due  to  a  combina¬ 
tion  of  systematic  errors  (errors  in  catalogued  star  positions, 
instrumentation,  and/or  adopted  longitude),  were  ignored. 
The  final  computation  of  ?7grav  used  station  elevations  inter¬ 
polated  from  the  ERTM2160  data  set1  (Hirt  et  al.  2014), 
and  included  the  ERTM2160  component  of  tj  to  reduce  the 
EGM2008  omission  errors.  However,  neglecting  the  effect 
of  elevation,  and  the  effect  of  ERTM2160  on  r\,  does  not 
appreciably  change  the  overall  statistics,  as  only  a  few  sta¬ 
tions  greatly  exceed  the  average  elevation  of  approximately 
200  m,  or  experience  highly  variable  terrain. 

For  all  BIH  stations  investigated,  the  longitude  differences 
A  -  Agrav  were  found  to  be  within  ±3  a  of  the  EGM2008 
commission  error  estimate  (Table  1).  The  weighted  aver¬ 
age  of  A  -  Agmy  was  statistically  no  different  than  zero 
(0.06"  =b  0.14"  1  a),  the  weights  being  the  reciprocal  val¬ 
ues  of  the  EGM2008-implied  commission  error  variances  of 
r\ gray  sec  cp.  Both  results  support  the  current  placement  of  zero 
longitude  as  being  due  to  the  local  deflection  of  the  vertical 
at  Greenwich,  with  no  change  in  orientation  occurring  at  the 
time  of  changeover  to  the  BTS  84. 

5  Other  conjectured  explanations 

Since  the  1960s,  realizations  of  zero  longitude  from  space- 
geodetic  reference  systems  have  tended  to  be  east  of  the 
Airy  Transit  Circle  by  an  amount  that  is  close  to  the  ver¬ 
tical  deflection  at  Greenwich  (Fig.  3).  This  offset  is  mostly 
a  consequence  of  artificial- satellite  operations  requiring  the 
plane  for  the  origin  for  UT1  to  pass  through  the  center  of  mass 
of  the  Earth.  In  1969,  Gebel  and  Matthews  (1971)  observed  a 
longitude  difference  of  5.64"  between  the  Airy  Transit  Circle 
and  zero  longitude  of  the  TRANSIT  navigation  system;  at  the 
time,  it  was  conjectured  that  DoV  at  Greenwich  contributed 
to  the  offset.  Relative  to  the  BIH  origin  for  UT1  ( e.g .,  seg¬ 
ment  OC "  in  Fig.  2),  TRANSIT’S  zero  longitude  was  further 
assumed  to  be  “quite  close,  i.e.,  about  1"  apart,  based  on  cer¬ 


1  Available  from  http://ddfe.curtin.edu.au/models/ERTM2160/. 


tain  reasonable  assumptions”,  although  the  details  behind  this 
assessment  were  not  published.  This  proximity  was  simply 
fortuitous:  the  original  longitude  reference  for  TRANSIT’S 
datum  was  implicitly  defined  by  a  regional  NAD  27  sur¬ 
vey  of  the  TRANSIT  antenna  site  (APLMND)  at  Applied 
Physics  Laboratory  (APL),  in  Laurel,  Maryland,  USA,  where 
the  longitude  difference  between  NAD  27  and  WGS  84  in  the 
vicinity  of  the  APL  is  about  1.2"  (Fig.  3),  as  estimated  by  the 
authors.  Later,  Dillon  et  al.  (1977)  adjusted  the  TRANSIT 
station  network  to  improve  compatibility  with  WGS  72. 

The  authors  have  also  identified  one  internet  article  (Mautz 
2004)  suggesting  DoV  as  a  probable  source  for  the  offset 
observed  with  GPS  from  Airy’s  meridian  line.  Nevertheless, 
the  early  research  of  Gebel  and  Matthews  (1971)  has  fueled 
an  increasingly  popular  opinion  that  APL’s  choice  of  longi¬ 
tude  reference  for  TRANSIT  promulgated  a  6"  offset  into 
subsequent  global  satellite  datums,  such  as  WGS  84  (IERS 
Reference  Meridian  2014).  In  reality,  the  early  geodetic  ref¬ 
erence  frames  of  TRANSIT  never  influenced  the  longitudinal 
orientation  of  WGS  84;  rather,  WGS  84  was  consistently  ori¬ 
ented  by  the  BIH  directional  origin  for  UT1  via  BTS  84,  a 
global  datum  independent  of  TRANSIT  (BIH  1985,  p.  B- 
6,  Table  3;  Defense  Mapping  Agency  1987;  Malys  1988). 
The  latest  realizations  of  WGS  84  are  aligned  with  the  ITRF 
established  by  the  IERS  (Wong  et  al.  2012;  NGIA  2014), 
which  also  evolved  from  BTS  84. 

Another  conjecture  for  the  longitude  offset  is  that  signifi¬ 
cant  errors  crept  into  the  BIH  timekeeping  process  gradually 
during  the  twentieth  century.  When  the  Airy  Transit  Cir¬ 
cle  began  to  experience  noticeable  problems  due  to  its  age, 
the  British  Time  Service  transitioned  its  timekeeping  pro¬ 
gram  to  other,  smaller  transits  after  July,  1927  (Gething 
1954).  The  Airy  Transit  Circle  was  already  out  of  regular 
timekeeping  service  when  the  BIH  began  using  its  “mean 
observatory”  to  determine  UT,  so  the  lack  of  direct  contri¬ 
bution  has  been  theorized  to  have  allowed  zero  longitude  to 
drift  away  from  Airy’s  transit.  However,  the  supplementary 
instruments  at  Greenwich  were  subject  to  the  same  vertical 
deflection,  and  continued  to  be  assigned  a  conventional  longi¬ 
tude  of  exactly  zero  by  the  BIH  (BIH  1933,  and  subsequent), 
with  their  observations  being  adjusted  for  their  longitude 
separation  from  Airy’s  instrument.  When  the  Royal  Green¬ 
wich  Observatory  (RGO)  astronomical  station  was  moved 
from  Greenwich  to  Herstmonceux  in  1957,  its  photographic 
zenith  tube  was  assigned  whatever  conventional  value  of 
longitude  gave  the  best  continuity  with  the  established  time 
services.  This  would  have  eliminated  the  discontinuity  nor¬ 
mally  expected  with  a  change  of  instrument  (Atkinson  1957). 
Although  some  of  the  processes  used  by  the  BIH  were  not 
transparent  to  users,  any  systematic  shifts  in  the  zero  longi¬ 
tude  were  apparently  limited  to  a  few  milliseconds  in  time 
at  most  (Mueller  1969;  Guinot  2000).  This  level  of  reported 
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Fig.  3  Selected  offsets  among  prime  meridians  from  the  “Space  Age”  and  the  Airy  Transit  Circle.  The  yellow  circle  provides  a  reference  for  each 
documented  offset 


accuracy  is  in  contrast  to  the  current  102  m  offset,  which 
would  be  more  than  l/$  s  in  time. 

Polar  motion,  which  has  been  well  measured  since  the  late 
nineteenth  century,  cannot  account  for  the  shift  of  the  prime 


meridian.  The  amplitude  of  the  periodic  part  of  polar  motion 
is  about  0.3  arc- second,  and  the  drift  of  the  mean  rotational 
pole  is  currently  about  0.5  arc-second/century  toward  the 
western  hemisphere  (IERS  2014).  The  effect  of  polar  motion 
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on  the  longitude  at  Greenwich  is  less  than  one  arc-second. 
Crustal  motion  from  plate  tectonics  is  also  much  too  small  to 
explain  the  observed  longitude  difference.  Britain  and  most 
of  Europe  are  moving  toward  the  northeast  at  about  3  cm/year 
(a  longitude  change  of  0.1  arc-second/century)  with  respect 
to  the  WGS  84  terrestrial  coordinate  frame  (Altamimi  et  al. 
2012).  The  accumulated  effect  relative  to  other  monitoring 
stations  amounts  to  only  a  few  meters  since  its  designation 
as  prime  meridian  in  1884. 

Improvements  in  the  fundamental  celestial  reference 
frame  and  in  astronomical  computations  for  positional 
astronomy  were  introduced  in  1984  after  10  years  of  plan¬ 
ning  (Kaplan  1981  and  references  therein),  but  there  is  no 
evidence  that  these  changes  affected  the  geodetic  longitude 
system  now  in  use.  These  changes  were  generally  of  the  order 
of  one  arc-second  at  the  changeover  date.  Although  small, 
the  changes  had  the  potential  to  affect  UT1  determinations 
and,  therefore,  the  longitude  system  of  BTS  84.  However,  the 
International  Astronomical  Union  specified  that  there  be  no 
change  in  the  value  or  rate  of  UT1  due  to  the  new  scheme 
(IAU  1977,  1980),  so  that  the  conventional  expression  for 
sidereal  time  in  terms  of  UT1  was  modified  to  provide  con¬ 
tinuity  in  the  UT1  determinations  (Aoki  et  al.  1982). 

6  Conclusion 

Zero  longitude  is  the  terrestrial  origin  for  Universal  Time 
(UT1),  a  realization  of  mean  solar  time  at  Greenwich  used 
to  define  the  rotation  angle  of  the  Earth  in  space.  The 
102-m  offset  between  the  Airy  Transit  Circle  and  zero 
longitude  indicated  by  a  GNSS  receiver  is  attributable  to 
the  fact  that  continuity  in  the  UT1  time  series  was  main¬ 
tained  in  the  BIH  reference  frames,  as  geodetic  longitudes 
replaced  astronomical  longitudes  when  space-geodetic  meth¬ 
ods  were  introduced.  This  continuity  condition  constrained 
the  BTS  84,  and  consequentially,  the  ITRF  zero  merid¬ 
ian  plane,  to  be  practically  parallel  to  the  orientation  of 
the  astronomical  prime  meridian  through  Airy’s  instrument 
that  is  aligned  to  local  gravity.  Extended  to  infinity,  these 
parallel  meridian  planes  sweep  past  the  same  stars  simulta¬ 
neously,  so  that  both  planes  indicate  the  same  astronomical 
time  (UT1).  The  difference  between  precise  GNSS  coordi¬ 
nates  and  astronomically  determined  coordinates  everywhere 
remains  a  localized  gravity  effect  due  to  the  deflection  of 
the  vertical;  thus,  no  systematic  rotation  of  global  longitudes 
occurred  between  the  former  astronomical  system  and  the 
current  geodetic  system. 

Because  the  value  of  the  DoV  at  Greenwich  and  the 
continuity  constraint  on  UT1  were  the  primary  factors  that 
influenced  the  current  location  of  the  ITRF  zero  meridian, 
the  1884  recommendation  of  a  prime  meridian  based  on  the 
“Observatory  in  Greenwich”  has  been  practically  maintained 


to  the  present  by  passing  the  zero  meridian  plane  through  the 
center  of  mass  of  the  Earth,  rather  than  through  the  center 
of  the  Airy  Transit  Circle.  Modern,  high-resolution  global 
gravitational  models  of  the  Earth  confirm  that  the  local  slope 
of  the  geoid  at  Greenwich  is  of  the  proper  sign  and  magni¬ 
tude  to  support  this  interpretation,  conclusive  to  within  the 
accuracy  of  the  deflection  model  (±0.5"  1  a  in  longitude  for 
the  EGM2008  model). 
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